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Abstract 

We examine packing of n congruent spheres in a cube when n is close but less 
than the number of spheres in a regular cubic close-packed (ccp) arrangement of 
[p^/2] spheres. For this family of packings, the previous best-known arrangements 
were usually derived from a ccp by omission of a certain number of spheres without 
changing the initial structure. In this paper, we show that better arrangements exist 
for all n ^ [p^/2] —2. We introduce an optimization method to reveal improvements 
of these packings, and present many new improvements for n ^ 4629. 


1 Introduction 

We consider the problem of finding the densest packings of congrnent, non-overlapping, 
spheres in a cnbe. Eqnivalently, we can search for an arrangement of points inside a nnit 
cube so that the minimum distance between any two points is as large as possible. The 
maximum separation distance of n points in [0,1]^ we denote by dn- 

To our knowledge, the optimality of dn is proved for n = 2, 3,4, 5, 6, 8,9 [TU], n = 10 
m and n = 14 [H]. Optimality of dn was conjectured for an inhnite family of packings 
where \p^/2] spheres are arranged in a cubic close-packed (ccp) structure |3]. We denote 
by g{p) = rP^/2] the number of spheres in these packings, with a maximum separation 
distance denoted by d^ = a/2/ (p — 1). 

In this paper, we examine arrangements when n is close, but less than g{p). For this 
family of packings dn = d^ is often assumed, to mean that the densest known arrangements 
are derived from ccp by omission of a certain number of spheres without changing the 
initial structure. Limiting values though were not provided. It was conjectured that dn is 
constant in the range 12 ^ n ^ 14 [3] and 29 ^ n ^ 32 [2]. A better packing was found 
for n = 12 [2]. Similarly, previous search results showed the same trend for 60 ^ n ^ 63, 
103 ^ n ^ 108, ..., 817 ^ n ^ 864 [3 13 Hg. 

In Section 2 we show that most of the listed packings can be improved by proving 
that dn > d'p for all n ^ g{p) — 2. We also provide a lower bound for these improvements. 
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In Section 3, we introduce an optimization method and improve the lower bound for 
4 ^ p ^ 21. We show that the described procedure can be used as a good packing 
method when n is slightly smaller than g{p). We run search to determine improvements 
for other packings when n = g{p) — r, for 3^r^6, r<p and 4 ^ p ^ 12. 


2 Existence of improved packings 


To simplify our notation, we will assume that the radius of all spheres in the packing is 
1 and that our task is to determine the smallest size of a cube that contains all spheres. 

We denote the family of all hnite sets of points such that the distance between any 
two points is at least 2 by 

-F = {S' C : ||si — S 2 II ^2 for all distinct Si, S 2 G S'} . 

Let Sn = {si,..., s„} G J-” and let Dc{si, sj) be the Chebyshev distance between any 
two points Si, Sj G Sn- The smallest edge length of a cube, with edges parallel to the axes, 
such that it contains all points S„ is equal to 

D{Sn) = max {Dc{si, Sj) : 1 ^ i < j ^ n} . 

We notice that the maximum separation distance dn can also be given as 

Theorem. The maximum separation distance of n points contained in a closed region 
bounded by a unit cube is larger than a/2/ (p — 1) for all n ^ |'p^/2"| — 2. 

Proof. Let Cp be a closed region bounded by a cube with an edge length Dp = 2/dp = 
(p — 1)a/2, where Cp is dehned by Cp = [0, Dpf. Let Cp be a set of g{p) sphere centers in 
a ccp arrangement, such that Cp E D, CpC. Cp and (0,0,0) G Cp (see Figure 1). 

For two sets of points A,BeD let 


h{A, B) = min {||a — h\\ — 2 : a E A,h E B} 

and let Lp = Cp\Cp-i (see Figure 2(a) for an example). We denote the improved packing 
of p(p) — 2 points by Pp = {si,..., Sg(p)_ 2 } G D such that 

Pp C Cp, (1) 

h{Pp, Lpj^i) > 0 . ( 2 ) 

Equations Q and (|^ directly imply that if these statements are true, then D{Pp) < Dp. 

Let Tp = T{Lp) be a set of points created by the translation of points from Lp by the 
function r such that 


Tp ED,TpC. Cp, h{Pp_i, Tp) ^ 0, 


( 3 ) 
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/i(Tp,Lp+i) > 0 (4) 

for all p > 2. If such a set Tp exists, then we can state that 

Pp = Pp-i U Tp, for all p > 2. 

We can easily show that P 2 exists by constructing it explicitly. If we want to give P 2 
in such a way to maximize the capability to translate points from L 3 and thus produce a 
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Figure 2: An illustration of arrangements L5, L5 1, L 5 2 and L5 3. 


better packing P 3 , we have to find the positive root of the polynomial 

+ 4a^ + Sa^ - 8 = 0, a > 0 ^ a = 0.818425 ... 

then P 2 is given as P 2 = {( 0 , 0 , a), ( 6 , 6, 0 )}, where b = — 0 ? j^. 

While Pp-i and Lp are separated ([^ , we can try to translate all points from Lp and 
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keep 0 , and Q true. We split Lp into three subsets, and perform the translations 
each of them while maintaining the given conditions. 

Let Lp = Lp^i U Lp ^2 U such that 

~ I ^ 2 , h) £ Lp : li, Ij ^ {p — 3)V2 for some distinct i,j G {1, 2, 3}| , 

Lp,2 = {(bj h, h) ^ Lp \ Lp i : ^3 > 0} , 

Lp,3 Lp \ Lp i \ Lp^2- 

For an illustration of sets L 5 , Ls^i, L ^^2 and Ls^s, see Figure 2. 

Now we can give Tp as 

Tp = Tp^i U Tp,2 U Tp^s, 

Tp,i = {{h - Ui^lTi{p),l2 -Ui^2Ti{p),h -Ui,3Ti{p)) : {h,l2j3) £ Lp^i} , 


on 




1 if Ij = {p — 1 ) \pl 
0 otherwise 


,M2j = 1,M3j = 


1 if Ij ^ 0 
0 otherwise ’ 


where 'ri(p), T2{p) and r^^p) are small numbers such that ri(p) > T2{p) > T3{p) and 
conditions (|^ and Q hold. If we additionally state that h(Pp_i,Tp i) = 0, h(Pp_i U 
Tp^i,Tp^ 2 ) = 0 and h(Pp_i U Tp i U Tp^ 2 ,^p, 3 ) = 0, we can give explicit solutions for ri, T 2 
and T3 as follows 

ri(3) = 272-2-0, 

n 2 (p) = ^ ^ + 272ri(p) + 4^ , for all p ^ 3, 

Mp) = ^ (^72r2(p) + 1 - \J-T2{py + 2V2 t 2 {p) + , for all p ^ 3, 

ri(p) = 72 + r 3 (p - 1) - \J-T3 {p - 1)^ + 2\/2 t 3 {p - 1) + 2, for all p > 3. 

We notice that D{Pp) = Dp — T^lp), hence the maximum separation distance can be 
given as 




(p- 1) 72 - r3(p) 
and while r 3 (p) > 0, then dn > 72/(p — 1). 

We denote the lower bound of improvements by 


, for all n ^ g{p) — 2 


□ 


Ip dp(p)_2 dp. 

By performing the calculations, we get particular values such as I 3 > 8.235 ■ 10“^^, 
I 4 > 1.276 ■ 10“^®,.... This is a rough approximation while condition Q needs to remain 
true for all p and does not allow us to further improve the packing for particular p. In 
the next section, we show that Ip is usually above this approximation. 

We did not hnd a way to improve the packings when n = g{p) — 1, and we conjecture 
that in this case dn = dp. 
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3 Optimization Approach 

Most of the existing packing methods focus on searching for a completely new arrangement 
of spheres, usually performing a search from a randomly given initial position of spheres. 
Such approach assumes that the packing of higher density can be reached after a certain 
number of iterations and multiple runs of the search procedure using different initial 
parameters The large number of iterations often limits the search procedure to 

the use of double or quadruple floating-point precision, to maintain computation speed. 
This precision is insufficient to detect improvements in many packings. Many of these 
approaches are adapted and modihed from widely known procedures for packing congruent 
circles in a square or a circle m la El El DS]. 

The method we suggest is based on a hypothesis that an improved packing can be 
reached just by the omission of two or more spheres from the ccp and by performing a 
translation of spheres using the available space made after we remove the spheres. 

The initial positions of the sphere centers we denote by Sp^r C Gp, as a set of g{p) —r 
points such that at least one sphere with a center s* G Sp^r can be continuously translated 
inside a cube container without overlapping with other spheres. More precisely, 

g{p)-r 

{o' G Cp : g ^ Sp^r-, {Sp^r \ {si}) U {g} G 7^ 0 . 

i=l 

We can see that the construction of Sp^r is possible for r ^ 2, and we experimentally 
determine solutions based on improvements reached for certain arrangements. To simplify 
the search procedure, instead of trying to hgure out the best performing arrangements 
Sp^r for each pair (p, r), we hnd removal patterns Rr = Gp\ Sp^r and use them to search 
for improvements for any p. 

Table 1 shows patterns in a simplihed notation where = | : s G We notice 

that only R 2 is the most likely an optimal pattern for all p > 1. 

r i?; 

2 {( 0 , 1 , 1 ), ( 1 , 1 , 0 )} 

3 {(0,0,0), (1,0,1), (2, 0,0)} 

4 {(0,1,1),(1,1,0),(1,1,2),(2,1,1)} 

5 {(0,1,1),(1,0,1),(1,1,0),(2,0,0),(2,1,1)} 

6 {( 0 , 0 , 0 ), ( 0 , 1 , 1 ), ( 1 , 0 , 1 ), ( 1 , 1 , 0 ), ( 2 , 0 , 0 ), ( 2 , 1 , 1 )} 

Table 1: Experimentally determined patterns i?) 

Using the initial arrangement Sp^r we try to perform the translation of each sphere 
using the limited set of translation vectors denoted by T. This algorithm can be described 
as follows: 
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. For a given initial set S'„ Sp^r repeat until D{Sn) < Dp-. 

. For each Sj G Sn do: 

. Randomly choose t E T, 

. Let V = {si + kt : k E [—1,1]}, 

. Let Vi = {qEv : {Sn \ {sj) U {q} E 7, D{Sn U {g}) = 

. Find endpoints a and b of the largest line segment ab C Vi such that s* G ab, 
. New position of Si is given as Sj <(— (a + b)/2. 

After we try to move all points from Sn, we say that we completed one iteration. 
Because of the very limited space to which we translate the spheres, and in order to 
minimize the number of required iterations, we usually set T = {0,1}^ \ {(0,0,0)}. We 
also tested performances when T takes different values such as { — 1,0,1}^, more or less 
reduced sets, but the improvements gained were always slightly worse. 

In practice, D{Sp^r) is slightly larger than Dp while coordinates are given with hnite 
precision. The described procedure stops when the hrst improvement is detected, but if 
we continue the search, we can improve the packing even more. It is also important to 
set the precision above the expected value of Dp — D{Sp^r), otherwise the improvement 
cannot be registered. Choosing the higher precision enables us to reach an improvement 
with less iterations, but only up to a certain level. We usually set the precision 1.5 times 
higher than the expected improvement. If the precision is too high, the search can be 
very slow, thus often we have to guess the range of possible improvements using a lower 
precision at hrst, and increase it if an improvement cannot be reached. 

This approach is different from the procedures used for sphere packing in the past, 
as in 121 El din], while we focus only on tiny changes/improvements in the high density 
structure. We cannot consider this approach as a good general packing method for r ^ p. 
Its main weakness is that, because of the small available space where spheres can be 
moved, random perturbations are hard to implement, or at least we did not hnd any 
good method to do it. Still, this method allows us to hnd improved packings in less than 
one second for some well examined cases even when high precision is not required, as for 
example n = 29, 59 or 60. 

The improvements attained are shown in the Tables 2 and 3. Table 2 lists values 
obtained for Ip with 4 ^ p ^ 21. Because of the slow computation times for p ^ 13, we 
ran a search with approximately 5000 iterations when improvement gains started to slow 
down. Table 3 lists other improved packings for 3^r^6,r<p and 4 ^ p ^ 12 using 
patterns described in Table 1. The results are listed as the best known values performed 
after a large number of iterations and multiple runs of the search procedure. 
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n 

P 

4 


30 

4 

7.34- 

10-68 

61 

5 

7.18 ■ 

10-80 

106 

6 

2.26 ■ 

10-314 

170 

7 

9.09 ■ 

10-622 

254 

8 

3.74- 

10-629 

363 

9 

7.51 ■ 

10-629 

498 

10 

5.00 ■ 

]^g-2584 

664 

11 

9.67- 

]^g-2563 

862 

12 

1.76 ■ 

]^g—4988 

1097 

13 

1.70 ■ 

]^g-5020 

1370 

14 

2.01 ■ 

]^g-5044 

1686 

15 

1.29 ■ 

]^g-5076 

2046 

16 

1.78 ■ 

]^g-5116 

2455 

17 

3.30 ■ 

]^g-5047 

2914 

18 

1.47- 

]^g-10118 

3428 

19 

4.20 ■ 

]^g-10121 

3998 

20 

1.16 ■ 

]^g-20344 

4629 

21 

2.46 ■ 

]^g-20582 


Table 2: Improved values of Ip 


n 

P 

r 

dn 



n 

P 

r 

dfi 

4 


29 

4 

3 

2.23 

■ 10- 

-12 


59 

5 

4 

1.95 

■ 10- 

-11 







60 

5 

3 

2.09 

■ 10- 

-20 







103 

6 

5 

3.38 

■ 10- 

-14 







104 

6 

4 

9.98 

■ 10- 

-47 







105 

6 

3 

1.34 

■ 10- 

-76 







166 

7 

6 

3.08 

■ 10- 

-21 

494 

10 

6 

4.77 

■ 10- 

-57 

167 

7 

5 

7.72 

■ 10- 

-31 

495 

10 

5 

6.96 

■ 10- 

-199 

168 

7 

4 

1.59 

■ 10- 

-87 

496 

10 

4 

1.72 

■ 10- 

-310 

169 

7 

3 

8.49 

■ 10- 

-148 

497 

10 

3 

1.28 

■ 10- 

-605 

250 

8 

6 

2.98 

■ 10- 

-28 

660 

11 

6 

2.83 

■ 10- 

-119 

251 

8 

5 

1.99 

■ 10- 

-43 

661 

11 

5 

1.99 

■ 10- 

-173 

252 

8 

4 

3.11 

■ 10- 

-102 

662 

11 

4 

1.40 

■ 10- 

-407 

253 

8 

3 

6.02 

■ 10- 

-153 

663 

11 

3 

3.47 

■ 10- 

-615 

359 

9 

6 

3.84 

■ 10- 

-28 

858 

12 

6 

8.40 

■ 10- 

-248 

360 

9 

5 

6.45 

■ 10- 

-44 

859 

12 

5 

2.68 

■ 10- 

-404 

361 

9 

4 

2.05 

■ 10- 

-101 

860 

12 

4 

7.74 

■ 10- 

-745 

362 

8 

3 

1.08 

■ 10- 

-152 

861 

12 

3 

2.91 

■ 10- 

-1212 


Table 3: Improvements reached using patterns Rr 
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